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o Effects of small-scale fluctuations in the neutrino radiation on core-collapse 

q supernova explosions are examined. Through a parameter study with a fixed 

+3 • radiation field of neutrinos, we find substantial differences between the results 

of globally anisotropic neutrino radiation and those with fluctuations. As the 
number of modes of fluctuations increases, the shock positions, entropy distribu- 
tions, and explosion energies approach those of spherical explosion. We conclude 
that global anisotropy of the neutrino radiation is the most effective mechanism 
of increasing the explosion energy when the total neutrino luminosity is given. 
This supports the previous statement on the explosion mechanism by Shimizu 
and coworkers. 

Subject headings: hydrodynamics — shock waves — starsmeutron — sup ernovae: general 



1. INTRODUCTION 

For many years since the first work of Colgate & White (1966), numerical simula- 
tions of core-collapse supernova explosions have been exciting topics. Until the beginning of 
the 1990s, almost all the simulations included the assumption of spherical symmetry (e.g., 
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Wilson 1985). Such one-dimensional simulations, however, were unable to explain the ob- 
served explosion energy and often failed to produce explosions (see, e.g., Liebendoerfer et 
al. (2001)). Inclusion of convective motion via the mixing length theory cures this prob- 
lem to some extent (Wilson & Mayle 1993; Bruenn, Mezzacappa, & Dineva 1995). However, 
spherical simulations based only on the Rayleigh- Taylor instability require considerably large 
initial fluctuations in density to explain the large-scale matter mixing. In addition, aspheri- 
cal explosion is also supported by the observation of SN1987A, where asymmetric ejecta are 
clearly observed (e.g., Wang et al. 2002). These lead us to multidimensional simulations of 
supernova explosions. 

At this time, the two- and three-dimensional simulations have been performed by several 
groups (Miller, Wilson, & Mayle 1993; Herant et al. 1994; Burrows, Hayes, & Fryxell 1995; 
Janka & Miiller 1996; Mezzacappa et al. 1998; Fryer & Heger 2000; Shimizu et al. 2001; 
Fryer & Wallen 2002; Kifonidis et al. 2003). In many multidimensional simulations, special 
attention is paid to the role of convection either near the surface of a nascent neutron star 
or in neutrino-heated regions above the neutrinosphere. It has been shown that large-scale 
mixing, caused by convection and convective overturn around the neutrino-heated region, 
increases the explosion energy and can trigger a successful explosion (Herant et al. 1994; 
Keil, Janka, & Miiller 1996; Janka & Miiller 1996). 

Because supernova progenitors such as OB stars are generally observed to be fast ro- 
tators (~ 200 km s^ 1 at the surface, P ~ 1 day; see, e.g., Tassoul 1978; Fukuda 1982), 
the resulting proto— neutron star can have a large amount of angular momentum after the 
gravitational collapse. Centrifugal force then deforms the rotating core into an oblate form. 
This will cause asymmetric neutrino radiation, in which the flux along the pole is enhanced 
over that on the equatorial plane. Janka & Monchmeyer (1989) first discussed the possi- 
bility of aspherical neutrino emission from a rapidly rotating inner core. They intended to 
evaluate the total neutrino energy outputs using the neutrino data of SN 1987A detected 
with the Kamiokande and Irvine-Michigan-Brookhaven experiments. It was argued in their 
paper that a neutrino flux along the pole might be up to a factor of 3 greater than that on 
the equatorial plane. 

Inspired by this work, Shimizu and coworkers (Shimizu, Yamada, & Sato 1994; Shimizu 
et al. 2001) proposed that the anisotropic neutrino radiation should play a critical role in 
the explosion mechanism itself, and carefully investigated the effects of anisotropic neutrino 
radiation on the explosion energy. They found that only a few percent enhancement in 
the neutrino emission along the pole is sufficient to increase the explosion energy by a large 
factor, and that this effect saturates around a certain degree of anisotropy. It should be noted 
here that the assumed rotational velocity of the inner core is very different between Janka & 
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Monchmeyer (1989) and Shimizu et al. (1994, 2001). Shimizu et al. (2001) concluded that 
the increase in the explosion energy due to anisotropic neutrino radiation occurs because 
cooling due to neutrino reemission is suppressed in anisotropic models. This is due to earlier 
shock revival and hence more efficient decrease of the matter temperature than those in 
spherical models. On the other hand, the neutrino heating itself is almost unchanged by the 
effect of anisotropic neutrino radiation. The neutrino heating dominates the cooling as a 
result, which increases the explosion energy, and leads to a successful explosion. 

In Shimizu et al. (2001), the geometric effects of neutrino radiation have been rigorously 
treated outside the neutrinosphere for the first time, although its flux on the neutrinosphere 
was assumed. Only a global form of anisotropy was assumed there; the maximum peak in the 
neutrino flux distribution was located at the pole, and the minimum at the equatorial plane. 
However, Burrows et al. (1995) have suggested that the neutrino flux can fluctuate with 
angle and time. Such fluctuations are due to gravitational oscillation on the surface of the 
proto— neutron star and have a completely different origin from that of globally anisotropic 
neutrino radiation. Thus, it is interesting to investigate how the small-scale fluctuations 
affect the explosion mechanism and compare the results with those of the global anisotropy. 

In this paper, we therefore introduce fluctuations in the neutrino flux in our numerical 
code by modifying the angular distribution of the neutrino flux. We aim to study the 
effects of these small-scale fluctuations on the shock position, the explosion energy, and the 
asymmetric explosion. Our numerical simulation is described in § 2, and the results are 
presented and discussed in § 3. Our conclusion is given in § 4. 



2. NUMERICAL SIMULATION 

Our simulation is performed by solving two-dimensional hydrodynamic equations in 
spherical coordinates. A generalized Roe's method is employed to solve the hydrodynamic 
equations with general equations of state (EOSs). The details of our numerical technique, 
together with the EOS and the initial condition used, are described in the previous article 
(Shimizu et al. 2001). In our study, we have improved the numerical code of Shimizu et 
al. (2001); the cells in the ^-direction were shifted by half of the cell size (Shimizu 1995) 
in order to avoid a numerical error near the pole, although the error was not serious for 
the investigation of the explosion energy. The computational region is divided into 500 
(r-direction) x 62 (^-direction) numerical cells. 

In the present paper, the local neutrino flux is assumed to be 

l u (r, 9) = \-JjT% (1 + c 2 cos 2 (M)) \, (1) 



-4- 



where a is the Boltzmann constant, and T v is the temperature on the neutrinosphere. In 
equation (1), the parameter c 2 represents the magnitude of anisotropy, and rig the number 
of waves in the ^-direction. The case of rig = 1 corresponds to the global anisotropy, namely, 
no fluctuation. We see in equation (1) that the neutrino fluxes in the x (equatorial) and 
z (polar) directions become l x = l„(r,9 = 90°) oc c\ and l z = l u (r,9 = 0°) oc ci(l + c 2 ), 
respectively. The degree of anisotropy l z /l x is then represented as 

Y = l + c 2 . (2) 



Note that equation (2) is different from that defined by Shimizu et al. (2001), (i z /i x )shimizui 
for ng — 1 and sufficiently small c 2 , we can relate them as c 2 ~ [(^/L)|himizu — 1] /2- 

The value of ci is calculated from given c 2 and rig so as to adjust the total neutrino 
flux to that in the spherical model. The total neutrino luminosity is obtained by integrating 
equation (1) over the whole solid angle, 

L u = J r%(r, 9)dQ = ^aT* 4tt Ci (l + c 2 §|^) , (3) 

which is equated to that of spherical explosion with the same T u , 

L7 = YQ aT * 47ri? NS- (4) 

In the above, i?NS is the radius of a proto— neutron star and fixed to be 50 km. By comparing 
equations (3) with (4), we obtain 

_ -"-NS /r\ 

1 l + c 2 (2ng-l)/(4n|-l)- K) 

It should be noted here that the magnitude of fluctuations in the neutrino flux dis- 
tribution for an observer far from the neutrinosphere (represented by c 2 ) and that on the 
neutrino-emitting surface (here we denote it as a) are different: the local neutrino flux is seen 
as equation (1) when we observe fluctuations on the surface of neutrino emission far from the 
neutrinosphere. It is preferable that we compare the results for the same value of a, since a is 
more directly related to explosion dynamics. The value of c 2 should, therefore, be calculated 
from a given a, depending on rig. Although it is difficult to calculate the exact relationship 
between c 2 and a, we can estimate it as follows. First, we assume that the strength of the 



- 5 - 



neutrino flux on the neutrinosphere is represented by a profile of step functions: 



for rig 



for rig 



for rig = 5 



/(*) 



f(0) 
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1 + a (1 > cos# > 1/2) 
1 (1/2 > cos# > 0), 

1+a (l>cos#>3/4, 1/2 > cos9 > 1/4) 
1 (3/4 > cos9 > 1/2, 1/4 > cos 9 > 0)) , 

1 + a (1 > cos 9 > 5/6, 2/3 > cos 9 > 1/2, 1/3 > cos 9 > 1/6) 
1 (5/6 > cos 9 > 2/3, 1/2 > cos 9 > 1/3, 1/6 > cos 9 > 0) . 



(6) 

We continue similarly for larger values of rig. The parameter a in equation (6) represents the 
magnitude of fluctuations in the neutrino flux on the neutrino emitting surface; the bright 
regions (/ ~ 1 + a) correspond to those where rising convective motion occurs, while the dark 
regions (/ ~ 1) correspond to those of sinking convection. We have assumed that the areas 
of bright and dark regions are the same. Note that the neutrino flux function is always bright 
(/ ~ 1 + a) when 9 = (along the pole) and dark (/ ~ 1) in the case of 9 = 7r/2 (on the 
equatorial plane), and that the fluctuations are essentially added to the global anisotropic 
model (ng = 1). 

The neutrino flux observed far from the neutrinosphere is obtained by averaging all con- 
tributions from the flux on the surface of the neutrinosphere, 1(0) = J d<p J d9f(9) sin 9 cos(9— 
O). Here is the inclination angle between the line of sight of an observer and the polar axis 
of the proto— neutron star. This means that fully geometric effects from an anisotropically 
radiating surface are included in neutrino radiation field. For example, there is a contribution 
from fluxes along the pole axis to those on the equatorial plane, which will tend to reduce 
an efficiency of anisotropy. The ratio of the local neutrino flux along the polar axis (l z ) to 
that on the equatorial plane (l x ) for an observer far from the neutrinosphere is described as 



/ d<j> \ f '(9) sin 6 cos OdO 
Jo Jo 

I d(p / /(#)sin#cos(#-7r/2)d0 
o Jo 

( 1 + 0-750a fn _ _ _ i 
1 + 0.391a tor Ue ~ X ' 



= < 



1 + 0.625a f — q 
1 + 0.438a torn 0-3, 

1 + 0.583a 



(7) 



I 1 + 0.457a 



for rig = 5, • 
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By comparing equation (7) with equation (2), we finally obtain the value of c 2 for each value 
of n e . 

In the following, we examine two model series: a = 0.31 (model series A) and a = 0.71 
(model series B). These values of a are chosen in such a way that the value of l z /l x for 
the global model (n e = 1) becomes 1.10 and 1.20, respectively. The values of c 2 for each 
fluctuation model (rig = 3,5) are accordingly calculated. These are summarized in Table 1. 
The neutrino temperature on the neutrino-emitting surface T v is assumed to be 4.65 and 
4.70 MeV. Note that the value of 1+a, 1.71, for the model series B roughly corresponds to 
the variation in the neutrino flux obtained by Burrows et al. (1995), which is at a factor of 
1.6. 



3. RESULTS AND DISCUSSION 

Figures 1 and 2 show the color-scale maps of the dimensionless entropy (Shimizu et al. 
2001) distribution with the velocity fields for the model of global anisotropy (rig = 1; models 
A1-T470 and B1-T470). At t = 82 ms after the shock stall, the shock front reaches r ~ 430 
km on the equatorial plane and r ~ 530 km at the pole for the model A1-T470. The shock 
front is prolate, since the neutrino heating along the pole is more intensive than that on the 
equatorial plane, resulting in a jetlike explosion. At a later stage (t = 244 ms), the shock 
wave is around a few thousand kilometers with large distortion. The entropy distribution of 
the model B1-T470 has a similar profile except that the shock front is more extended. 

In Figures 3 and 4, the results for models with fluctuations (models A3-T470, A5-T470, 
B3-T470, and B5-T470) are depicted. When rig = 3, we find that the degree of asymmetry 
is smaller than that of the global anisotropy. In particular, the shock position for these 
models is less extended than the globally anisotropic one. This trend becomes more obvious 
for the model of c 2 = +0.035 and rig = 5 (model A5-T470), where the shock front is almost 
spherical and its radius is only about 1300 km. In the case of the model B5-T470, the shock 
front is distorted because of a strong hydrodynamic flow along the pole, although this does 
not affect the explosion energy (see discussion on the energy later in this section). 

The profile of the explosion energy is found to be closer to that of spherical explosion 
as the mode number of fluctuations increases. Figure 5 shows the evolution of the explosion 
energy, as well as the thermal, kinetic, and gravitational energies for the models of T v = 4.70 
MeV. The difference between the globally anisotropic model and the models with fluctuations 
is prominent: the energy gain for the case of rig = 1 is the highest among others at all stages 
of the explosion. It is also seen that the explosion energy decreases as the mode number of 
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fluctuations in the neutrino flux increases and finally approaches that of spherical explosion. 
For series B, we have obtained similar results. We observe that asymmetry in explosive 
motion is more enhanced than that for each model of series A. 

In Figure 5, we compare the results of the explosion energy for the two model series. 
No significant difference is found between the two, although the result of the model Al- 
T470 becomes larger than that of the model B1-T470 at the later stages of the explosion. 
It has been shown (Shimizu et al. 2001) that the explosion energy increases as the degree 
of anisotropy becomes larger for not so large a degree of anisotropy, and finally saturates 
at ( l z I l x ) shimizu ~ 1-2. Our new result shows that the final explosion energy of the more 
anisotropic model B1-T470 is smaller than that of the model A1-T470 (see Fig. 5). The 
difference may be attributed to the fact that the assumed forms of the local neutrino flux 
are different (compare eq. [1] here with eq. [5] in Shimizu et al. 2001). In the present paper, 
we have assumed the form of the neutrino fluxes that has more sharply concentrated flux 
on the pole. Therefore, the neutrino heating and rising convection are focused on the pole, 
and those in the equatorial direction are extremely reduced. The shock wave of the model 
B1-T470 appears to be too weak on the equatorial plane at t ^ 300ms, which causes a energy 
loss (Shimizu et al. 2001). Such features are clearly seen in the entropy distribution of the 
model B1-T470. 

The effectiveness of global anisotropy becomes more pronounced as T v is decreased. 
Figure 6 shows the same energy evolution as Figure 5, except for T v = 4.65 MeV. The 
difference between the model with global anisotropy and those with fluctuations is extremely 
remarkable. The globally anisotropic models succeed, while all the other fluctuated and 
spherical models fail to explode except for the model B3-T465. Note that difference in the 
neutrino temperature is only 1%, which indicates a sensitivity of the supernova problem on 
the neutrino luminosity and energy. Note also that an increase of the explosion energies 
of the models with fluctuations at t ~ 500 ms is physically meaningless, because we have 
not taken into account a decay of the neutrino luminosity when t > 500 ms (e.g., Wilson & 
Mayle 1988; see also Shimizu et al. 2001). 

We found that there are remarkable differences in the explosion energy depending on 
the mode of the fluctuations and that larger number of modes in the fluctuations makes 
the result closer to that of spherical explosion, irrespective of the model series A or B. Any 
small-scale fluctuations on the neutrinosphere are greatly averaged out when the neutrino 
emission is observed far enough from the neutrino-emitting surface. Moreover, we found 
that a certain broad space is needed to be heated by neutrinos to revive the stalled shock 
wave rigorously and that the global anisotropy (n e = 1) is the most effective to increase the 
explosion energy. Burrows et al. (1995) suggested that the neutrino flux can fluctuate not 
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only with angle but with time. Such time fluctuations are expected to reduce further the 
efficiency of anisotropy, which needs to be confirmed in the future. 

4. CONCLUSION 

We have investigated the effects of small-scale fluctuations in the neutrino flux on the 
core-collapse supernova explosion. In order to examine the effect of the degree of anisotropy 
itself on the explosion, we have studied two model series parametrically. We specified the 
neutrino radiation field taking its geometric effects into account in each model. We found that 
the global anisotropy (ng = 1) and the local fluctuations (rig > 1) in the neutrino flux have 
quite different effects on the explosion mechanism, that is, the shock dynamics, the explosion 
energy, and the explosion asymmetry. Since the small-scale fluctuations are averaged out for 
radiative and hydrodynamic reasons, the results including fluctuations become closer to that 
of spherical explosion. Consequently, the global anisotropy is the most effective mechanism of 
increasing the explosion energy. Note here that the explosion energy could differ substantially 
depending on the neutrino temperature (the difference between 4.70 and 4.65 MeV is only 
1%). This indicates that the supernova problem is very sensitive to the neutrino energy 
and luminosity. However, the total luminosity cannot be simply increased to explain the 
observed explosion energy because such treatment leads to the problem of Ni overproduction, 
especially in the case of essentially spherical models. We therefore conclude that globally 
anisotropic neutrino radiation is of great importance in actual supernova explosions. This 
supports the claim made by Shimizu et al. (2001). 

The global anisotropy can originate from rotation of a proto— neutron star or a hot 
spot on the neutrino-emitting region, while the small-scale fluctuations are considered to be 
resulted from gravitational oscillation or uniform convection. It will be very interesting if 
any evidence of anisotropic neutrino radiation is observed at facilities like Super Kamiokande 
(Hirata et al. 1987; Suzuki 1998) and SNO (Poon et al. 2001), together with detailed optical 
observations (e.g., Wang et al. 2002). 

We are grateful to the anonymous referee for useful comments and suggestions that 
improved this paper. 
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Table 1. Simulated Models 



Model Series 


a a 


Model 


n e h 


c 2 c 


T/(MeV) 






A1-T465 


1 


0.100 


4.65 






A1-T470 


1 


0.100 


4.70 


A 


0.31 


A3-T465 


3 


0.051 


4.65 






A3-T470 


3 


0.051 


4.70 






A5-T465 


5 


0.035 


4.65 






A5-T470 


5 


0.035 


4.70 






B1-T465 


1 


0.200 


4.65 






B1-T470 


1 


0.200 


4.70 


B 


0.71 


B3-T465 


3 


0.101 


4.65 






B3-T470 


3 


0.101 


4.70 






B5-T465 


5 


0.068 


4.65 






B5-T470 


5 


0.068 


4.70 



a Magnitude of fluctuations on the neutrino emitting sur- 
face. 

b Mode number of fluctuations. 

c Magnitude of fluctuations far enough from the neutrino 
emitting surface. 

d Temperature on the neutrinosphere 
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colorscale & velocity field colorscale & velocity field 




x [kin | x [km J 

max 2.7x10° min 5.3xlO- E max 4.6x10° mm 2.0x10-' 

Fig. 1. — Color-scale maps of the dimensionless entropy distribution and the velocity fields 
for the model of n 9 = 1 and T v = 4.70 MeV of series A (model A1-T470). Left: t = 82ms 
after the shock stall, Right: t = 244ms. 




Fig. 2. — Same as Fig.l, except for series B (model B1-T470). Left: t = 82ms after the 
shock stall, Right: t = 249ms. 
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colorscale &: velocity field eolorscale & velocity field 




max 4.8x10° rnin 1.9x10-' max 4.0x10° mill 1.7X10 1 



Fig. 3. — Same as Fig.l, except for the case of fluctuated neutrino flux. Left: n e = 3 (model 
A3-T470) at t = 254ms after the shock stall, Right: n e = 5 (model A5-T470) at t = 250ms. 




Fig. 4. — Same as Fig. 3, except for series B. Left: ng = 3 (model B3-T470) at t = 257ms 
after the shock stall, Right: uq = 5 (model B5-T470) at t = 247ms. 
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Fig. 5. — Evolution of thermal and kinetic energy (E th + E kin ), gravitational energy (-E grav ) 
and explosion energy (E exp \) for the models of T v =4.70 MeV. Solid line corresponds to the 
case oi no = 1, short-dashed line rie = 3, long-dashed line no = 5, and dotted line C2 = 
(spherical). Left: model series A, Right: model series B. 




Fig. 6. — Same as Fig. 5, except for T v = 4.65 MeV. Left: model series A, Right: model 
series B. 



